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We study e + e~ — > J/ipcc by measuring the invariant mass distribution recoiling against fully 
reconstructed J/ip decays, using 124 fb -1 of data collected with a center-of-mass energy of 10.6 GeV 
with the BABAR detector. We observe signals for ri c (lS), XcO, and r] c (2S) in the recoil mass distribu- 
tion, thus confirming previous measurements. We measure a(e + e~ — + J/ip +cc) B(cc — > > 2 charged) 
to be 17.6±2.8(stat)+2;i(syst) fb, 10.3±2.5(stat)+j'*(syst) fb, and 16.4±3.7(stat)+3 *(syst) fb with 
cc = r] c (lS), XcO, and rj c (2S), respectively. 

PACS numbers: 13.66.Bc, 12.38. Qk, 12.38. Bx, 14.40.Gx 



Prompt J/ip and ip(2S) production in e + e~ annihila- 
tions around yfs — 10.6 GeV has been observed by both 
the BABAR p| and Belle 0] experiments. These interac- 
tions provide an opportunity to study both perturbative 
and non-perturbative effects in QCD and to search for 
new charmonium states 0,0. 

Belle @ reported the observation of r) e (lS), XcO, an d 
i] c (2S) in the mass distribution of the system recoiling 
against a reconstructed J/ip in e + e~ annihilations. The 
production cross sections measured by Belle are about 
one order of magnitude higher than those predicted by 
non-relativistic QCD (NRQCD) calculations 0,110 for 
e + e~ — > 7* — > J/ip cc reactions, where cc is a char- 
monium state with even C-parity. There have been 
attempts @ H 0, 0, E3 to reconcile the large dis- 
crepancy between the observed cross section and pre- 
dictions, and the validity of NRQCD approximations 
has been questioned [3, 113] • It has also been suggested 
that at least part of the double charmonium production 
might be due to two virtual-photon interactions [lOj, i.e., 
e+e~ — > 7*7* — > J/ip cc, where odd C-parity states could 
be produced. Belle updated its observation and explored 
the origin of the J/ip cc events • 

In this paper we present a measurement of the cross 
sections for e + e~ — * J/ip rj c (lS), e + e~ — ► J/ip XcO, an d 
e+e~ — > J/ip rj c {2S), and set limits on the yields for other 
known charmonium states produced in association with a 
J/ip . We calculate the mass (M roc ) of the system recoiling 
against a fully reconstructed J/ip via: 

M^={^s-E*j^f-pf^, (1) 

where y/s is the e + e~ annihilation energy in the center- 
of-mass (CM) system, and Ej^ and p*j^ are the energy 
and momentum of the J/ip candidate in the CM system. 

In this paper, we analyze 112 fb -1 of data collected at 
the peak of the Y(AS) resonance and 12 fb -1 at y/s — 
10.54 GeV, just below the T(45), with the BABAR de- 
tector jlj operating at the asymmetric energy PEP-II 
e + e~ storage ring. The BABAR detector includes a five- 
layer, double-sided silicon vertex tracker (SVT) and a 
40-layer drift chamber (DCH) in a 1.5-T solenoidal mag- 
netic field, which detects charged particles and measures 
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their momenta and specific ionizations (dE/dx). Photons 
and electrons are detected with a CsI(Tl)-crystal elec- 
tromagnetic calorimeter (EMC). An internally reflecting 
ring-imaging Cherenkov (DIRC) is used for particle iden- 
tification. Penetrating muons are identified by an array 
of resistive-plate chambers (RPC) embedded in the steel 
of the flux return (IFR). 

We select events with at least five well reconstructed 
charged tracks in the DCB, within the fiducial volume 
0.41< 8 <2.54, where 6 is the polar angle. Electron can- 
didates have a pattern of specific ionization (dE/dx) in 
the DCH, a Cherenkov cone angle, an EMC shower en- 
ergy divided by momentum, and a number of EMC crys- 
tals that are consistent with an electron hypothesis. A 
muon candidate is selected on the basis of energy de- 
posited in the EMC, the number and distribution of hits 
in the IFR, and the match between the IFR hits and the 
extrapolation of the DCH track into the IFR. A more de- 
tailed explanation of particle identification is given else- 
where [jj . 

A pair of oppositely charged lepton candidates origi- 
nating from a common vertex is selected as a J/ip can- 
didate if its mass (m(£ + £~)) falls within [-50,30] MeV/c 2 
(for e+e") or [-30,30] MeV/c 2 (for £i + /-0> of the nomi- 
nal J/ip mass of 3.097 MeV/c 2 [Hi- In the calculation of 
m(e + e~), electron candidates are combined with nearby 
photon candidates in order to recover some of the energy 
lost through bremsstrahlung radiation. These mass in- 
tervals are referred to as the J/ip mass windows. In order 
to improve the P*j^ resolution, we perform a kinematic 
fit where the J/ip candidate is constrained to have the 
nominal J/ip mass. 

There are two main background sources in this anal- 
ysis: events with genuine J/ip mesons and combinato- 
rial background. The region 60MeV/c 2 < \M (£+£-) - 
M(J/ip)\ <200MeV/c 2 , defined as the J/ip mass side- 
bands, where M(J/ip) is the nominal J/ip mass, is used 
to estimate the combinatorial background due to random 
tracks. This background is largely rejected by particle 
identification, and by a requirement on the lepton hclic- 
ity angle in the J/ip decay, |cos6>;| < 0.9, as shown in 
Fig- □» and He). 

The largest backgrounds are due to real J/ip mesons 
from QED processes such as J/ip or ip(2S) mesons pro- 
duced via initial state radiation (ISR). J/ip mesons from 
B meson decay have p* < 2 GeV/ c and do not consti- 
tute a background for recoil masses below 6.6GeV/c 2 . 
Most QED backgrounds have low multiplicity, and may 
have electrons or photons escaping detection along the 
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beam line. These backgrounds are suppressed by the re- 
quirement of at least five charged tracks and the follow- 
ing requirement: for each event we calculate the energy 
deposited in the EMC plus the energy that can be at- 
tributed to an undetected electron or photon, 



E, 



QED 



E 



EMC ± Pn 



(2) 



where -Eemc is the total energy deposited in the EMC, 
and p m i ss is the missing momentum in the lab frame in 
the event. We require Eqed — Eb cams < — l.OGeV as 
shown in Fig. ^b) and ^d) , where -Ebcams is the sum 
of the e + e~ beam energies calculated in the lab frame. 
We reject the J/ip background from ip(2S) events by ve- 
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FIG. 1: Distributions of (a) cosf?; and (b) Sqed — Sbeams 
in the data, (c) cos#; and (d) -Eqed — Bbeams in the signal 
Monte Carlo. The arrows point to where the selection criteria 
are applied. 

toing events if the invariant mass of the J/ip candidates 
combined with any pair of oppositely charged tracks with 
pion mass hypothesis is within 15MeV/c 2 of the ip(2S) 
mass. 

The recoil mass distribution for events in the J/ip mass 
window is shown as points with error bars in Fig. [21 The 
ISR ip(2S) background is estimated using a Monte Carlo 
sample of ISR ip{2S) events. The ip(2S) feeddown back- 
ground from continuum production is estimated using 
continuum ip(2S) events selected in the data. 

The spectrum in Fig. [2] is fit to the sum of signal 
functions representing the rj c (lS), XcQ, and rj c (2S) line- 
shapes, plus a second-order polynomial background func- 
tion. The signal line shapes are obtained by convolut- 
ing the Breit-Wigner line shape of each resonance with 
a fixed-width Gaussian representing the recoil mass res- 
olution function. The widths of the Gaussians are deter- 
mined from a Monte Carlo simulation of the momentum 
of the reconstructed J/ip ; the J/tp momentum resolution 
is different for the J/tp — > e + e~ and J/ip — ► sam- 
ples, but independent of the recoiling system. This shape 
in turn is convolved with a long radiative tail that is cal- 
culated to 0(a 2 ) 01 f° r ISR photons that carry off an 
energy greater than 10 MeV. The free parameters in the 
data fit are the coefficients for the background parame- 
terization, the event yields for each resonance, the masses 
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FIG. 2: The fit to the recoil mass distribution is represented 
by the solid curve. The dashed curve is a second-order polyno- 
mial representing the background. The points with error bars 
refer to the events in the J/ip mass window. The histograms 
represent different sources of backgrounds. 



TABLE I: Result of the fits to the recoil-mass spectrum. The 
errors are statistical only. Where indicated, the value of the 
corres pon ding parameter is fixed to the current world av- 
erage [16J. The primary fit is obtained including signals of 
T) C (1S), XcO, and r/ c (2S). The event yield for the other reso- 
nances is determined by including each resonance in the pri- 
mary fit. 



Recoil 
System 


Number 
of Events 


Mass 
(MeV/c 2 ) 


Total Width 
( MeV/c 2 ) 


Vc(lS) 


126 ± 20 


2984.8 ± 4.0 


fixed 


XcO 


81 ± 20 


3420.5 ± 4.8 


fixed 


Vo(2S) 


121 ± 27 


3645.0 ± 5.5 


22 ± 14 


J/4* 


-26 ±13 


fixed 


fixed 


Xol 


-5 ±16 


fixed 


fixed 


Xc2 


-12 ± 16 


fixed 


fixed 


4,(2S) 


30 ±27 


fixed 


fixed 



of the resonances, and the rj c {2S) total width. The to- 
tal widths for the r] r (lS) and the \c0 are fixed to their 
world average values [l(| of 17.3 MeV/c 2 and 10.1 MeV/ c 2 , 
respectively. The fit is performed simultaneously to the 



recoil mass spectra in the J/ip — > e + e~ and J/ip — > 
samples, and the total event yield for each resonance is 
given by the sum of the yields in each mode. 

The fit result is given in Table ^ and is shown as the 
solid curve in Fig. [21 Other known charmonium states 
may also be produced in association with the J/ip via 
two virtual-photon interactions. We therefore attempt 
to include in our primary fit each one of the other known 
charmonium resonances in turn to determine their event 
yields, which are presented in Table [I] We find no evi- 
dence for J/ip, Xcij Xc2, or ip{2S) in the mass spectrum 
of the system recoiling against a J/ip . 

The topological branching fraction is unknown for the 
ri c (\S), Xco, and rj c (2S), so we report the product of the 
branching fraction for final states with more than two 
charged tracks (2?>2(cc — ► > 2 charged)) times the dou- 
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TABLE II: Summary of systematic errors: variations of cross 
sections and masses due to the selection and fitting procedure 
(Fit), particle identification (PID) efficiency, and recoil-mass 
scale uncertainty. A M refers to the mass difference between 
the r) c (2S) and 7] C (1S). 





Variations(%) 


Variations( MeV/c 2 ) 




in Cross-section 




in Mass 


Source 


Vc(lS) 


XcO 


Vc(2S) 




XcO 


r) c (2S) AM 


Selection 


+3.5 
-8.3 


+0.3 
-9.2 


+ 12.6 
-15.6 


+3.0 
-0.2 


+ 1.2 

-0.7 


+ 1.0 7 9 
-5.1 '' z 


Fit 


+6.7 
-8.1 


+13.5 
-14.2 


+6.8 
-8.3 


+0.1 
-3.4 


+9.9 
-8.0 


+3.3 +7.1 
-3.7 -2.3 


PID 


±3.5 


±3.5 


±3.5 








Mass Scale 








±1.5 


±1.5 


±1.5 


Sum 


+8 
-12 


+ 14 
-17 


+ 15 
-18 


+4.5 
-5.0 


+ 11.5 

-9.5 


+4.9 +7.1 
-7.8 -7.6 



ble charmonium production cross section. In order to 
include the effect of ISR, the yields reported in Table |U 
are calculated with a line shape based on a model of the 
t/s dependence of double charmonium production model. 
To allow a direct comparison of experimental results, we 
follow the same method used by Belle 0] to remove this 
model dependence by determining cross section values 
that correspond to the non-tail fraction of the fit shape 
(frad = 0.61) where no ISR photon with an energy 
greater than 10 MeV is radiated. We use 



J /if) cc) B 



>2 



Nccfrad 



B(J/ip ^i+l-)Ce ' 



(3) 



where N c5 is the event yield, C is the integrated luminos- 
ity, B(J/ip — > £ + £~) is the J/ip branching fraction, and e 
is the detection efficiency The value of e is determined 
using a Monte Carlo simulation with the assumption that 
exclusive J/ipr) c (nS) production is P wave and that ex- 
clusive J/ip XcO production is S wave, as expected for 
a single virtual-photon process. The efficiency is deter- 
mined to be (28.8±0.7)% for the rj c (lS), (31.5±0.7)% for 
the XcO, and (28.9±0.8)% for the T] C {2S). 

The systematic error is estimated taking into account 
contributions from the event selection and the fitting 
procedure, the particle identification efficiency, and the 
recoil-mass scale uncertainty. The contributions from un- 
certainties in integrated luminosity and J/ip branching 
fraction are negligible. The contributions from individ- 
ual sources (listed in Tabic ITT|) are added in quadrature, 
except for the systematic errors due to the mass-scale 
uncertainty, which are added linearly, to determine the 
total systematic errors. 

We obtain a(e + e~ — > J/ip cc)B(cc — > > 2 charged) to 
be 17.6 ± 2.8±2;? fb for J/ip r) c (lS), 10.3 ± 2.5t\i * for 
J/ip XcO, and 16.4±3.7i|o fb for J/ip r? c (2S). Throughout 
this paper, the first error is statistical and the second sys- 
tematic. Our values of the cross sections are consistent 
with Belle's measurements ^3 f° r all three resonances. 
The cross sections measured by both experiments are 
much larger than those predicted by many NRQCD cal- 



TABLE III: Comparison of cross-sections (a x £?>2 in fb) with 
Belle's results [l4|l . and with theoretical expectations that do 
not include the £>>2 factor. 



J/ip cc 




XcO 


Vc(2S) 


BABAR 


17.6±2.8±2:i 


10.3 ± 2.5±i;| 


16.4±3.7lg;o 


Belle [14] 


25.6 ±2.8 ±3.4 


6.4 ± 1.7 ± 1.0 


16.5 ± 3.0 ± 2.4 


NRQCD [6J 


2.31 ± 1.09 


2.28 ± 1.03 


0.96 ±0.45 


NRQCD [4] 


5.5 


6.9 


3.7 



culations. 

From the fit to the recoil mass spectrum we determine 
the r) c (2S) mass to be 3645.0 ± 5. 5±|;| MeV/c 2 , and the 
total width to be 22 ± 14 MeV/ c 2 . The systematic errors 
are mainly due to the uncertainty on the J/ip momen- 
tum measurement. We use ISR J/ip and ISR ip(2S) data 
samples to determine the momentum shifts away from 
the expectations for ISR events. Assuming a constant 
momentum shift, we obtain the recoil mass uncertainty 
for J/ip cc processes due to the J/ip momentum uncer- 
tainty. The mass difference (AM) between the rj c (2S) 
and r] c (lS) does not significantly depend on the absolute 
momentum scale and common systematic errors mostly 
cancel. We measure AM = 660.2±6.8+7;gMeV/c 2 , which 
is in good agreement with the mass difference previously 
reported by this experiment [Tsj and by other experi- 
ments [illfl9| . 

In summary, we have measured the cross sec- 
tion for double charmonium production cr(e + e~ — > 
J/ip cc) B(cc ^> 2 charged) for J/ipr] c (lS), J/ip XcO, 
and J/tpr] c (2S). We confirm the unexpectedly large 
cross sections previously reported by the Belle exper- 
iment for these processes. No evidence is found for 
e+e" -> J/ipJ/ip, J/ipXd, J/^Xc2, or J/^^(2S). We 
also measure the mass difference between the r/ c (2S) and 
the T) C (1S). 
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